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We study the dependence of the spin splitting on the number TV of atomic layers, using first- 
principles calculation for Au(lll) surface. When the slab of the atomic layers is sufficiently 
thick, the lower split state has a minimum away from f , which is known as the Rashba 
effect. As the number of layers decreases, the minimum approaches f , and it is located at 
f for N < 14. This crossover is analyzed in detail using two models: a tight-binding model 
and a bilayer nearly-free-electron model. It is demonstrated that the features of surface band 
dispersion are clearly understood as the result of the competition between the interference 
of the surface states on both sides and the spin-orbit interaction. 
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1. Introduction 

Recent progress in spintronics has been increasingly evoking interests in the role of spin- 
orbit interaction in condensed-matter physics. Spin-orbit interaction induces many intriguing 
phenomena, such as the spin-Hall effect, topological insulator and magnet oelectric effect. At 
surfaces, it causes spin splittings even in nonmagnetic materials. This effect was pointed out 
by Rashba. 1 ) He analyzed a two-dimensional electron with a momentum k moving freely on 
the xy plane under an electric field along the z direction: 

k 2 

H R( k ) = — + aa-(e z x fc), (1) 
2m 

where a are the Pauli matrices. The second term originates from spin-orbit interaction. The 
constant a in the isotropic Rashba Hamiltonian, eq. (1), is proportional to the gradient of the 
potential along the z direction: a = 4m 2 c 2 ^- The dispersion relation of this Hamiltonian has 
two branches for the eigenvalues 

E ± (k) = f-±ak. (2) 
2m 

The two eigenstates for nonzero k have opposite spin polarization. The spin splitting is ex- 
plained by the lack of inversion symmetry. Kramers' theorem states that in a crystal invariant 
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under time reversal operation —io y K^ where K is the complex conjugation operator, two 
Bloch states having opposite wave vectors and opposite spins are degenerate: E^ — E_^i- 
Independently of this, a system having inversion symmetry satisfies E^ a — E_^ a for a =t, I- 
Hence in a crystal with both time reversal and inversion symmetry, E^ = £74 must be satis- 
fled, that is, spin splitting is prohibited. On a surface, inversion symmetry is broken and spin 
splitting can occur. 

Experimental studies on the Rashba effects reported so far are not only for sur- 
faces of metals [e.g., Au(lll), 2 ' 3 ) Ag(lll), 3 ) and Bi(lll) 4 )] but also for surface alloys 
[e.g., Ag/Au(lll), 5 " 7 ) Bi/Ag(lll), 8 ) and Sb/Ag(lll), 9 )] including the binary surface alloy 
Bi a: Pbi_ x Ag(lll) 10 ) and semiconductor heterostructures, 11 ' 12 ) focusing particularly on how 
the magnitude of the splitting could be tuned by varying the element and formula. Among 
these, the Au(lll) surface is a well studied system. 

LaShell et al 2 ^ observed Rashba splitting, 110 meV at the Fermi level, of the L-gap surface 
states on Au(lll) via an angle resolved photoemission spectroscopy (ARPES) experiment, 
which is much greater than the zero- field spin splitting for a two-dimensional electron gas, 
which is on the order of meV in semiconductor heterostructures. 11 ' 13 ) Their results have been 
theoretically confirmed using a tight-binding model 14 ) and first-principles electronic structure 
calculations. 3 ' 15 " 17 ^ The anisotropic Rashba splitting on Au(110) has also been studied by 
Simon et a/. 18 ) with the k p perturbation theory by introducing an anisotropic Rashba term 
H^ ls (k) = \ x k x a x + \yk y (jy. Petersen and Hedegard 14 ) pointed out that, in a monolayer tight- 
binding model for Au(lll) (where the system has spatial inversion symmetry), the constant of 
the Rashba Hamiltonian, eq. (1), must include the effect not only of the inherent intraatomic 
spin-orbit interaction but also of the transfer integrals 7 = (p z (R)\V\p n (R + x))(n — x,y) 
between p orbitals deformed by the surface potential for the reproduction of the splitting. 
Nagano et al 17 ^ demonstrated more quantitatively with first-principles calculations that the 
spin-orbit interaction occurs only in the vicinity of each atom 19 ) and that the asymmetric shape 
of the squared wave function ^(z)] 2 along the normal of the surface is crucial for the large 
spin splitting. The asymmetry comes from the mixing of atomic orbitals of different parity, 
and the energy splitting is proportional to the integral of the potential gradient multiplied by 
\i/j(z)\ 2 , of surface states on Au(lll). It is in contrast with Ag(lll) and Sb(lll), for which 
|^(z)| 2 's are less asymmetric. The importance of asymmetric feature of surface wave functions 
was also pointed out for Bi/BaTiO 3 (001). 20 ) 

So far, the Rashba effect has been studied for solid surfaces. However, because of the 
recent progress in nanoscale fabrication technology, thin films on substrates are becoming 
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target materials of researches on systems where spin-orbit interactions play important roles. 
The possibility of tuning the thickness of an overlayer on a substrate induces interests in 
the thickness dependence of the Rashba splitting 5-7 ) from a practical viewpoint since the 
controllability of the Rashba splitting is desirable for the development of spintronic devices. In 
the present work, we investigate the magnitude of the spin splitting as a function of the number 
of layers using both first-principles and model calculations. The present work also provides 
technical insights for an electronic structure calculation using a finite slab. A sufficiently thick 
slab is quite often used, mimicking a semi-infinite surface system in an ordinary electronic 
structure calculation, because it can be treated as a three-dimensional periodic system when 
it is separated from its duplicated image by vacuum regions. Electronic structure calculations 
for semi-infinite systems without using a slab reported so far are few due to the complicated 
formalism and the technically difficult implementation of such calculations. 21 ' 22 ) We find that, 
when the slab is thick enough, the energy dispersion of the surface state is well described by 
eq. (2). The eigenvalue of the lower branch decreases as the wave vector moves away from T. 
We find that, as the slab gets thinner, the energy gap at f becomes larger, and a crossover 
takes place. The energy splitting is the Rashba splitting when the slab is thick, whereas it is 
characterized by hybridization between the surface states of both sides of the slab when it is 
thin. 

2. Computational Details 

First-principles electronic structure calculation is based on the density functional theory 
(DFT). We adopt the projector augmented- wave (PAW) method 23 ) using the Quantum MA- 
terials Simulator (QMAS) package 24 ) within the local-density approximation (LDA). 25 ) We 
implemented fully relativistic calculation using two-component pseudo Bloch wave functions 
on QMAS. The total energy of the system is calculated as a functional of the 2x2 density 
matrix defined as 

occ. 

Paa'ir) = ^2ipnk*(r)*ip nka >(r), (3) 

n,k 

where a, a' = a, /3 are spin indices. It is in contrast to spin-independent nonrelativistic or 
scalar relativistic DFT calculations, in which charge density is the fundamental quantity. The 
formulation of fully relativistic electronic structure calculations using two-component wave 
functions has been proposed for norm-conserving pseudopotentials, 26 ) ultrasoft pseudopoten- 
tials, 27 ' 28 ) and PAW formalisms 29 ) in detail. In a fully relativistic calculation, noncollinear 
magnetism and spin-orbit interaction can be naturally introduced. 
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We solve the Dirac equation as follows. Under a central electrostatic potential, a four- 
component energy eigenfunction solution of the Dirac equation is of the form 

*^r) = l( f ^'^). (4) 

i(jjj Z K is not an eigenstate of L 2 , but a simultaneous eigenstate of the operators H, J 2 ,J Z ,S 2 , 
and /C, whose eigenvalues are E,j(j + l),j z ,3/4, and ft, respectively. J = L + S is the 
total angular momentum operator. For the definition and properties of /C in detail, readers 
are referred to, e.g., Ref. 30. ft can take only j + 1/2 or —(j + 1/2). For ft = +{j + 1/2), 
/a = j =t 1/2 and Zb = j T 1/2. is the spinor spherical harmonics 31 ) defined as 



l±jz±l/2 




for j = Z ± 1/2, which is a simultaneous eigenstate of the operators J 2 , J Z1 L 2 , and 5 2 . /j« and 
(/j^ are the solutions of the following system of differential equations: 32 ) 

^(E-V + 2c 2 )g jK = 0, 

[i-^ + - c {E-V)f jK = 0. (6) 

In ordinary condensed matter physics, the upper two components are much larger than the 
lower two components in eq. (4). We therefore drop the lower components, and ^jj zK becomes 
two-component and now an eigenfunction of L 2 . Thus the elimination of gj K from eq. (6) 
leads to the single differential equation to be solved for the construction of a fully relativistic 
potential: 

^d 2 1(1 + 1) 1 dV ( d n 



2M(r)(E-V) 



fin = 0, (7) 



dr 2 r 2 2M{r)(? dr I dr r 
where we denote I a by / and M(r) — 1 + (E — V)/2c 2 . The term involving ft originates 
in the spin-orbit interaction and causes j-splitting of the energy spectrum. Since eq. (7) is 
independent of j z , it has solutions degenerate in 2j+l = 21+2 and 21 for j = 1+1/2(k = —l — l) 
and j — I — 1/2(k — /), respectively, for fixed /. j-averaged ft is hence 

( K ) = (-/ - 1) • + 1 • = -1. (8) 

w v ; 4/ + 2 4/ + 2 v ; 

If we replace ft with (ft) in eq. (7), the scalar relativistic equation is obtained, whose energy 
eigenvalues depend only on /. By introducing ft(A) = —1 — Xl and — 1+A(Z+1) for j = 1+1/2 and 
j = I — 1/2 for common /, respectively, we can continuously move from the scalar relativistic 
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(A = 0) equation to the fully relativistic (A = 1) equation by varying the strength A of the 
spin-orbit interaction. 

In the present study pseudo wave functions are expanded in plane waves with an energy 
cutoff of 30 Ry for 6 x 6 fc-points for the surface Brillouin zone (see Fig. 1(a)), which give 
sufficiently converged results for the purpose of the present study. The slab is constructed by 
stacking Au atomic layers. The lattice constant is fixed to the experimental value 4.078 A of 
bulk Au, 33 ^ and the slab is separated from its duplicated image in the neighboring unit cell 
by a vacuum region of 15 A width. The atomic positions are not relaxed. 

3. Results and Discussion 

3.1 First-principles calculation 

Figure 1(b) shows the fully relativistic electronic band structure of a 22-layer slab. There is 
a bunch of states in the energy region from —8 to —1 eV. They have a strong 5d character. The 
bands having 6s and 6p orbital characters cross the 5d bands. Hybridization between them 
leads to the sd-derived surface states near —7.5 eV and the sp-derived surface states 16 ' 17 ) near 
—0.5 eV at f . The latter correspond to the experimentally observed 2 ) L-gap surface states 
with the spin splitting. The experimental and calculated Fermi wave vectors for the inner and 
outer L-gap surface states are shown in Table I. Our values are in good agreement with the 
experimental values. 

Table I. Experimental and calculated Fermi wave vectors (in A -1 ) for the inner and outer Z/-gap 
surface states. 









Lout 


Exp. 


LaShell et a0 


0.153 


0.176 




Reinert et al u "> 


0.167 


0.192 




Nicolay et al 3 ^ 


0.172 


0.197 


Calc. 


Henk et a/. 15 ) 


0.149 


0.172 




Mazzarello et al 16 ^ 


0.159 


0.191 




Present Work 


0.172 


0.201 



We then changed the number of layers and carefully examined the features of the band 
structures of the L-gap surface states in the vicinity of f . Figure 2(a) shows the band dispersion 
of the L-gap surface states at V for N = 13, 16, and 19. There are two branches, each of which 
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r m k r 



Fig. 1. (Color online) (a) Surface Brillouin zone of Au(lll) and primitive reciprocal lattice vectors, 
(b) Fully relativistic electronic band structure of a 22-layer slab. The origin of energy is set to the 
Fermi level, (c) is a blowup of (b), showing the vicinity of the spin-split L-gap surface states at T. 

is twofold-degenerate. (There are four states in total: two degrees of freedom from spin, and 
the other two from the number of surfaces.) The calculated energy band of the inner surface 
states was found to be a monotonically increasing curve as a functions of wave vector for all 
the TV's. It was found, however, that for N = 16 and 18, the band takes a minimum value at 
a nonzero wave vector fco, while for N = 13 the outer surface state band has a minimum at F. 

In order to look into this behavior, we artificially changed the strength A of the spin-orbit 
interaction and observed the variation of the band structure. The result for N = 19 is shown 
in Fig. 2(b). We can see that the bands of the outer surface states reach their minimum only 
at r for A = and 0.5, while the minimum is observed off T for A = 0.75 and 1. This clearly 
indicates that the spin-orbit interaction is crucial for this behavior. 

We performed systematic calculations for the number of atomic layers N = 11-22, and 
found that the minimum is away from T for N > 15. Figure 2(c) shows a plot of ko and the 
energy gaps AE at both f and fco as functions of N. It is seen that AE(T) rapidly decreases as 
N is increased, while the variation of AE(ko) is much smaller despite the outward movement 
of fco away from f . 
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In the following two subsections, we analyze these results in more detail using model 
calculations. It is elucidated from a bilayer model that the difference in the features of these 
bands, which has the critical number of layers, 15, comes from the competition of the strengths 
of the spin-orbit interaction and the interference between the surface states on both surfaces. 

3.2 Tight-binding model 

Before proceeding to the analysis of a bilayer model, we examine here the interference of p- 
derived surface states on the uppermost and lowermost layers via a tractable spin-independent 
tight-binding calculation by varying the number of layers of a slab. 

Within the tight-binding formalism, the Hamiltonian matrix for a wave vector k — k\b\ + 
k 2 b 2 ly m g on the surface Brillouin zone (see Fig. 1(a)) is given by 

Hij(k) = ^2e ikR t i0j R: (9) 

R 

where the sum runs over the lattice points corresponding to the surface unit cell and tiojR 
is the transfer integral between the orbital i in the home unit cell and the orbital j in the 
unit cell at R. The orbital indices represent p x ,P y , and p z orbitals in the present study. We 
take into account only the transfers between nearest-neighboring atoms, t p and t Sl which are 
often conventionally denoted by (ppn) and (ppcr), respectively. The explicit expressions for 
the intralayer Hamiltonian matrix in bulk, Hq, are thus easily calculated: 

3 1 

Ho xx (k) = -(c 2 + c 12 )t p + -(4ci + c 2 + c 12 )t 8 

1 3 

H 0yy(k) = -(4ci + C 2 + C 12 )t p + -(c 2 + C 12 )t s 

H 0zz (k) = 2(ci + c 2 + c 12 )t p 
"\/3 

Ho X y(k) = H 0yx (k) = ~y(-C2 + c 12 )(-t p + t s ) 

H 0yz (k) = H 0zy (k) = H 0zx (k) = H 0xz (k) = 0, (10) 

where q = 2 cos 2nki(i = 1, 2) and c\ 2 = 2 cos 2ir{ki+k 2 ). Those for the interlayer Hamiltonian 
matrix in bulk, Hi, are also calculated: 

Hi xx {k) = ^(3 + 3ei + 4e 12 )tp + ^(1 + ei)t s 
H iyy( k ) = + lle i + 8e u )t p + -^(1 + e x + 4e i2 )t 5 



1 

Vs, 

12 



H Uz (k) = -(1 + ei + e 12 )(tp + 2t s 



Hi xy (k) = H lyx (k) = —(-1 + ei)(t p - t s 
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a/2 

H lyz (k) = H lzy (k) = — (1 + ei - 2e 12 )(t p - t a 

\/6 

Hi zx (k) H lxz (k) — (1 - ei)(t p - t s ), 



(11) 



where ei = e _z27rfcl and e±2 = e _ ^ 27r ( fcl+/c2 ). We have set the p orbital energies to zero. 

We incorporate the surface perturbation on the orbital energies by adding a constant 
£, which is described by the surface effect Hamiltonian Hse = diag(£, e, e). The 37V x 3N 
Hamiltonian matrix for a TV-layer slab is thus given by 

(h q + h se h x \ 



H(k) 



Hq Hi 

H\ Ho H x 



V 



(12) 



Hq Hi 
H\ Hq + HseJ 

The band dispersion of the slab is obtained by solving the time-independent Schrodinger 
equation H{k)c^ — Ej^Ck, where c& is the 3iV-dimensional column eigenvector and is its 
corresponding energy eigenvalue. The bulk band structure is obtained as a continuum of real 
E such that real 9 exists satisfying det(if|e~^ + H - E + Hie 19 ) = 0. 

We do not incorporate in Hse the transfer integrals 7 between p orbitals deformed by the 
surface potential, which were introduced by Petersen and Hedegard 14 ) for a monolayer system. 
The reason for this is as follows. We study multilayer systems. The translational symmetry is 
inevitably broken at the surface and thus the perturbed transfer integrals are not necessary 
for surface states. Furthermore, the matrix elements coming from 7 vanish at f and hence the 
following discussion would be unchanged even if they were incorporated. 

Let us inquire into the surface states at f . We rearrange the order of the bases as 
{\pl), . . . , |pf ), \pl), . . . , \p*), \pl), . . . , |pf )}, where \p\) is the Bloch sum of the p { orbital 
on the Z-th layer for the wave vector k. The Hamiltonian matrix of the slab at f then looks 
block-diagonal: 



H(T) 



H J 



(13) 
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l 



\ 



(14) 



for i =||, -L and the transfer parameters as tjj = 3t p + 3t S: tQ = 6t p ,t\ = 5t p /2 + t s /2, and 
= t p + 2£ s . Two iifll in H(T) correspond to the twofold energy spectra. One is ^-derived 
and the other is p^-derived. gives the p z -derived spectrum. H l is of the same form as 
the tight-binding Hamiltonian matrix for a finite one-dimensional chain only with nearest- 
neighbor transfers between N sites. Such a system has been analyzed in detail by Davison and 
Grindlay. 35 ) This one-dimensional system admits two surface states when \z\ > l(z = s/t\). 
The eigenvectors corresponding to the surface states are, in the limit of N — >> oo, arbitrary 
linear combinations of degenerate localized states given by 

CL„ = A JV (sgnz) n e -^- 1 ) 

c Rn = A N (sgnzTe^ n - N ^ 



(15) 
(16) 



H = In |^|, A N = e^ N -^ 2 . 



sinh n 



(17) 



sinh jiN 

cl and cr are the solutions decaying from the first and N-th sites into the bulk, respectively, 
which have the common eigenvalue E = £q + (sgnz)2t\ cosh/x. For finite iV, however, the 
bonding and antibonding states are formed by the surface states localized on the individual 
surfaces, leading to an energy gap AE. As an example, Figs. 3(a) and (b) show the band 
structure of a 15-layer slab with t p = —0.3 and e — —2.5 in units of t s . Two purely p z -deiived 
surface states appear below the bulk band continuum at f , whose wave functions are localized 
on both surfaces, as shown in Fig. 3(c). 

Let us estimate AE by constructing the normalized bonding and antibonding states as 



(18) 



V2(1±S)' 

We calculate the gap as a difference in the expectation values of the Hamiltonian between 
them and obtain 

2 



AE = t c-H i c- - t c + H i c + 



(19) 



1-S 2 

Since t c\ J H t CYi is expected to contain contributions mainly from several sites around the 
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midpoint of the chain, this term is on the order of e~^ N_1 \ while ES is obviously on the 
order of Ne'^-^. Equation (19) can thus be reduced to 

AE^2ENe~^ N ~ 1 \ (20) 

which is consistent with the rapid decay of the gap at T obtained in the first-principles 
calculation (see Fig. 2(c)). 



3.3 Bilayer model 



We here propose a bilayer two-dimensional free-electron model that continuously connects 
the two extreme situations, in which either the interference of the surface states or the spin- 
orbit interaction is so strong that the other is negligible. This model consists only of two layers 
to each of which free electrons are confined, feeling the spin-orbit interaction and the interfer- 
ence from the other layer. The individual layers are described by the Rashba Hamiltonian, eq. 
(1). The Hilbert space for this model is spanned by four states, namely, free-electron-like wave 
functions on the upper and lower layers with up and down spins: {|U f), |U j), |L t) 3 |L i)}. 
The Hamiltonian matrix in k space reads 



/ 



H(k) 



JsL 

2m 



-a(ik x + k y 



a(ik x - k y ) 
V k e~ i<t>k 




2m 



V k e %<t>k 






V k e i<t>k 



\ 



id 

2m 



-a(ik x - k y ) 



a(ik x + k y 

2m 



(21) 



/ 





where m is the effective mass of an electron confined to the layer. The component V^e 1 ^ 
(V& and 0& are real) is responsible for the interference between the surface states on the 
confronting layers. This matrix has two eigenvalues, to be compared with eq. (2), 



a 



2 k 2 , 



(22) 



each of which is twofold degenerate regardless of the k dependence of the interference compo- 
nent. Considering the basic knowledge of the matrix theory that eigenvectors belonging to the 
same eigenvalue can be freely rotated by an arbitrary unitary matrix, we adopt the following 
orthonormalized eigenvectors of H{k)\ 




k±L) 



\ ) 




(23) 
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where d& = a/Vk measures the strength of the spin-orbit interaction compared with the 
interference between the layers. |VVh(-)u) and |^fe+(-)L) belong to E + (_^(k). The degener- 
ate eigenvectors for the same branch have opposite expectation values of the spin operator 



(S)k±\J = —(S)k±L = (Tkydki ±k x dk, — l)/2y 1 + & 2 d|, which is natural because of the in- 
version symmetry of the system. The ratios of the weights on the layers of the eigenvectors 
are given by 

^^ = §^ = 1 + 2^, (24) 
where P\j and Pl are the projection operators onto the upper and lower layers, respectively. 
Equation (24) implies that, for d& ^ 0, \ifjk±u) an d \^k±h) have major components on the 
upper and lower layers, respectively. Let us verify that dy. continuously connects the two 
limits, the strong interference limit and the strong spin-orbit interaction limit. In the former 
limit, that is —> 0, the eigenvectors are |^fe±u) = t (0, 1, 0, ±e~^ fc )/v / 2 and \^k±h) = 
H±e^,0,l,0)A/2. |^ +( _)u) and |^ +( _) L ) are antibonding (bonding) states with equal 
weights on both layers with purely down and up spins, respectively (see Fig. 4(a)). In the 
latter limit, on the other hand, that is d& — » oc, the eigenvectors are |^fed=u) — t ( z pi(kx — 
ik y )/k, 1, 0, 0)/y/2 and \ipk±h) = t (0, 0, 1, =p(^z + ik y )/k)/y/2, for which the two layers are 
decoupled and each of these layers becomes an ordinary two-dimensional free-electron system 
with the Rashba term (see Fig. 4(b)). 

We then examine the dispersion relation of the bilayer system, regarding the interfer- 
ence component to be real and fc-independent, Vb, for easier qualitative understanding. The 
eigenvalues, eq. (22), take extrema at k = when Vb ^ 0, giving an ^-independent energy 
gap of E+(0) — E-(0) = 2|Vb|. In addition, if |Vb| < ma 2 = V C1 E_ takes a minimum at 
ko = \/V 2 — Vq jot and thereby the gap is E+(ko) — E_(ko) = 2V C , independent of Vq (see 
Fig. 4(c)). The Vb-independence of this gap is consistent with the result of the first-principles 
calculation, shown in Fig. 2(c). The ratio K/|Vo| characterizes the crossover of the band split- 
ting from the bonding-dominant nature to the spin-dominant nature. Combining the result 
of eq. (20) for the tight-binding model analysis, we can obtain the condition for the number 
of layers in order for the spin-split bands to be reproduced in a slab calculation. The gap at 
k = should be identified with AE = 2ci7Vc^~ 7V , where c\ and C2 is independent of N, and 
hence the condition for nonzero ko is 

Cl Nc^ N < V c . (25) 
This is a necessary condition for reproducing the spin-split band structure characteristic to 
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the ideal Rashba system, described by the Rashba Hamiltonian, eq. (1). It tells us explicitly 
that a slab consisting of sufficiently many layers allows the Rashba splitting to be reproduced 
despite the presence of the inversion symmetry in a practical electronic structure calculation, 
that is, the smallest N that satisfies this condition is 15 in our first-principles calculation for 
a Au(lll) slab. 

4. Summary 

We performed fully relativistic first-principles calculations of the electronic structure of 
Au(lll) slabs with varying the number N of atomic layers and the strength A of spin-orbit 
interaction. The variation of the features of the bands of the sp-derived surface states in the 
vicinity of f was examined with either one of N and A varied and the other one fixed. We 
found, for N > 15, that the finite wave vector ko exists, for which the energy gap is nearly 
unchanged when N is increased, in contrast to the gap at f . The rapid decay of the gap at f 
as a function of N was explained by the simple tight-binding calculation of a slab including 
only p orbitals. We adopted the bilayer two-dimensional free-electron model and were able 
to clearly understand the qualitative behavior of the first-principles band structure with N 
and A varied. In addition, combining the results of the tight-binding and bilayer models, we 
obtained the explicit condition of N and a for occurrence of spin-split bands in an electronic 
structure calculation using a slab. From the success of the bilayer model, the combination of 
first-principles calculations and simplified low-dimensional models is expected to help one to 
understant real systems more complicated than a slab, such as an overlayer on a substrate and 
an interface. Furthermore, the insights obtained in the present work are useful for achieving a 
reliable electronic structure calculation of a surface system because they help one to distinguish 
artifacts coming from the finiteness of a slab and the intrinsic properties of the surface system. 
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Fig. 2. (Color online) (a) Fully relativistic band structures in the vicinity of the L-gap surface states 
at f for TV = 13, 16, and 19. The origins of energy are set to the respective Fermi levels, (b) shows 
the band structure for N = 19 with varying the strength A of the spin-orbit interaction, (c) plots 
ko and energy gaps at f and ko as functions of N. ko is the finite wave vector, for which the outer 
surface states take a minimum band energy. 



13/16 



J. Phys. Soc. Jpn. 



Full Paper 




Fig. 3. (Color online) (a) Tight-binding band structure of a 15-layer slab. The shaded region is the 
bulk band continuum. Energies are in units of t s . (b) is a blowup of (a), showing the vicinity of 
two purely p^-derived surface states at f . The wave functions and weights of these surface states 
are shown in (c). 



(a) v k ^w (b) v k ^0 




Fig. 4. (a) and (b) are schematic illustrations of a bilayer two-dimensional free-electron model in 
Vk oo and Vk limits, respectively, drawing the wave functions of the four eigenstates. 
Filled arrows represent spin polarization, (c) shows band structures of the bilayer system with Vb 
varied. 
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